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We show how to construct the measure of the path integral in lattice gauge theory. This measure 
contains a factor beyond the standard Haar measure. Such factor becomes relevant for the calcula¬ 
tion of a single transition amplitude (in contrast to the calculation of ratios of amplitudes). Single 
amplitudes are required for computation of the partition function and the free energy. For U(l) 
lattice gauge theory, we present a numerical simulation of the transition amplitude comparing the 
path integral with the evolution in terms of the Hamiltonian, showing good agreement. 
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I. INTRODUCTION 


In lattice gauge theory it is customary to compute the expectation value of an observable 0\U], such as 
the Wilson loop, which is given by a ratio of Euclidean path integrals 


< Ul,t = -|-oo|0[[/]|fl, t = —oo >= 


J[dU] d[U]] exp[-S[U]/h] 
J[dU] exp[-S[U]]/h] 


( 1 ) 


Here [dU] denotes the Haar measure of the group of gauge symmetry, like for example, C/(l), SU{2), SU{3), 
and U denotes the link variables, being elements of such group. The state O denotes the vacuum. The 
reason why it is customary to compute such ratio of path integrals is due to the Monte Carlo method with 
importance sampling (like Metropolis) which works only for such ratio. In contrast, let us consider a single 
transition amplitude, like 

< Uf,,t = T\U,n,t = 0 >= J[dU] exp[-5[t/]/fi]|^f;J . (2) 

The states \Uin >, \Ufi > denote a Bargman link states, which are defined by assigning a value Uij to each 
link ij on the lattice in a fixed time slice. Now the measure [dU] is no longer given by the Haar measure 
only, but there is a factor involved for a lattice of iV — 1 intermediate time slices. It is the objective of 
this article to construct such measure. Physical scenarios which require the calculation of single amplitudes 
are: amplitudes of decay reactions, the partition function the free energy F{j3) given in terms of the 

partition function, or matrix elements involved in the construction of the Monte Carlo Hamiltonian P, |2| . 


II. MEASURE IN 1-D QUANTUM MECHANICS 

Consider a Lagrangian of the form 

L(x,x) = — V{x) . 


( 3 ) 
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The transition amplitude in imaginary time, expressed as a path integral over TV — 1 intermediate time 
slices is given by 


(cc/il exp[-iVT/ft]| 
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y^in — ^0: ^fi — ^Ny ^0 — T/N . 


(4) 


The measure 


d/r = dxi 


dxN-i 



(5) 


contains the physical parameters of mass to, Planck’s constant h and the length of a time slice ag, but is 
independent of the potential. It can be obtained by computing the propagator of the kinetic term for a 
single time slice 


exp[- 


P^ao 

2mh 


\Xr. 



m 

2agh 


{Xfi 



( 6 ) 


This propagator has the following properties. First, when ag goes to zero, the propagator goes over to 
6{xfi — Xin)- The r.h.s. actually is a representation of the <5 function by a Gaussian. Second, for any value 
of ag and Xm one has 




(7) 


Its interpretation is that the free propagator is the solution of a stochastic process (diffusion), and represents 
the probability for a random walker starting at Xi„ at t = 0 to arrive at x/i at time t = ag. Summed over 
all possible hnal destinations, the probability must be one. 


III. MEASURE IN LATTICE GAUGE THEORY 


Let us consider QED on the lattice without fermions. The groim of gauge symmetry is U{1). We keep 
in mind how to construct a lattice Hamiltonian (see e.g, CreutzQ, Rothej^) via the transfer matrix, by 
splitting the lattice action into terms of time-like and space-like plaquettes, 

S[U] = 4- E [l-7?e(C/n)] + 4- E [ 1 " ^^(C/n)] 

g^ag ^ a ^ 

•—•time —Hfce ^space — like 

= S''’’“[C/]-f S'P°*[C/] . (8) 


Like in quantum mechanics, the extra factor in the measure is determined solely by the kinetic term. Thus 
we consider the single transition amplitude 

< C//,|exp[-iL'=“T/fi]|I/„ > , (9) 

where denotes the kinetic term of the Kogut-Susskind lattice Hamiltonian Q, given by 
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( 10 ) 
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The Hamiltonian requires to choose a gauge and the Kogut-Susskind Hamiltonian has been obtained using 
the temporal gauge (Uume-Uke = !)■ The propagator, Eq.@, expressed as path integral reads 


(C//,|exp[-i?'=“T/;i]|C/, 


= / Z[dU] exp 
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( 11 ) 


The amplitude factorizes into independent amplitudes for each spatial link ij. In order to compute the 
factor Z, let us consider a single link and its time evolution for a single time step (T = ao). 


(H/,|exp[-H-'=*’^ao/;i]|C/„) = Z exp 


dg^ao 


[1 COs(^Cyfi 


( 12 ) 


Like in quantum mechanics also in lattice gauge theory the Euclidean propagator has a probabilistic 
interpretation, and the analogue of Eq. 0 holds. 


dC//,(H/,|exp[-Lr'=“oo/?i]|C/,„) = 


da 




27r 


Z exp 


hg'^ao 


[1 - cos(a/i - a^n)\ 


= 1 


Defining A = a/{hg^a()), and using the Bessel function of imaginary argument 

1 r 

Io{z) = — d0exp[zcos(0)] , 

Jo 

Eo. ltT!^ yields 

Z(A) - 

loiA) ■ 

Eor fixed T, lim^v^cxj means Oq ^ 0 and A —> oo. Then one finds 

Z{A) = V^[l - in-i - + 0(H-3)] . 


(13) 


(14) 


(15) 


(16) 


In the limit oq ^ 0 the leading term Z{A) = V27rH is sufficient to guarantee that the amplitude in Eq. 113 
goes over to S{Ufi — Uin), as should be. In quantum mechanics the leading term Z{A) = ^Aj (27r) with 
A = ml{haQ) is the exact result (see Eq.@). However, in lattice gauge theory, the sub-leading terms are 
important and can not be neglected. This has been confirmed by a numerical simulation discussed below. 


IV. COMPARISON OF PROPAGATOR FROM PATH INTEGRAL AND HAMILTONIAN 
TIME EVOLUTION - A NUMERICAL SIMULATION 

Let us consider a lattice which in spatial direction has a single link ij and in time direction has N time 
slices. We compare the propagator expressed via Hamiltonian time evolution with the path integral, using 
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the measure of the group integral taking into account the previously calculted factor Z. 

(t//,|exp[-iJ'=*'^r/ft]|17„) 
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The propagator expressed by Hamiltonian time evolution is given by 


(17) 



*^final ^initial 


FIG. 1: Transition amplitude for single link state as function of initial and final link, Uin = exp(iQi„) and Ufi = 
exp(iQj'i). Transition time T = 1, lattice spacing a = 1. Comparison of Hamiltonian time evolution, Eg. 11811 (bold 
line) with path integral, En. ll 71 (thin lines, number of time-slices N = 8, 16, 32,..., 256). 


< H/il exp[-i7''“r/S]|t7i„ >= ^ exp[-^^^u^] cos[n(ai„ - a/i)] . (18) 

This has been calculated using the basis of eigen states of the electric field operator 

iij\Xij >= Xij |Aij > , Xij = 0, ±1, ±2,... (19) 

and the connection from the Bargmann link basis to the electric field string basis given by the scalar 
product 


< A|17 >= (t/)^ . (20) 

The path integral, Eq. 03 , has been evaluated by numerical integration using 10000 mesh points in each 
time slice. The results as function of the angle o/i — cxm and of the number of time slices N is shown 
in Fig.^. One observes convergence when increasing N = 8, 16,32,..., 256. The relative error, shown in 
Fig.El varies between 10 ^ and 10 ^ depending on the angle and decreases monotonically when increasing 
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FIG. 2: Relative error of transition amplitude shown as function of = 2, 4,8,..., 256 and oifi — cein- 


N = 2,4, 8,, 256. On the other hand, when taking into account only the leading term in Ea. (ll6ll the 
relative error was found to be in the order of 10 — 20%. 
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